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1. Introduction 


Holographic duality has provided condensed matter physicists with a novel way to study the properties 
of strongly correlated transport and dynamics in theories without quasiparticles [1, 2, 3]. Charged black 
holes become dual to states of matter at finite charge density and temperature. At finite density, however, 
it becomes crucial to include the effects of momentum relaxation in order to obtain a finite direct electrical 


1 







conductivity, cJdcj at zero frequency and momentum.^ These momentum relaxing effects may either be 
studied by numerically constructing a charged black hole which breaks translational invariance [4, 5, 6], or 
analytically, so long as the black hole breaks translation invariance only weakly [7, 8, 9, 10], translational 
invariance is broken via massive gravity [11, 12, 13, 14], “helical lattices” [15, 16] or “Q-lattices” [17, 18, 19]. 

In this paper we show, through a direct holographic computation, that the low-frequency conductivity 
is given by a Drude peak in a wide class of holographic metals where translational symmetry is weakly 
broken: 


a{uj) 


(^dc 

1 — \U)T 


( 1 ) 


These metals have coherent transport without quasiparticles, in the language of [20]. Furthermore, fide 
and r are identical to the results computed - independently of holography - through the memory function 
formalism [21, 22, 23]. Memory functions have been used, in conjunction with holographic methods, to 
evaluate conductivities in the past [7, 24], as well as in non-holographic condensed matter models for 
strange metals as well [25, 26, 27, 28]. More recently it has been noted that exact computations of 
(Tdc “ the direct current conductivity at zero momentum and frequency - admit temperature scaling in 
agreement with results found through memory functions [8, 9]. In fact, the correspondence hinted at in 
these works is exact, including all other prefactors. A short computation within the memory function 
framework necessitates the emergence of a Drude peak (though with complications when magnetic helds 
and/or strong charge diffusion is present [29]). It is not obvious from a purely holographic computation 
why this Drude peak should be universal, and how it should arise - other than resorting to the argument 
that these metals should be describable by the memory function formalism. We will address these points 
in this paper. 

We begin with a review of why the Drude peak is universal, via hydrodynamic arguments, in Section 
2. We then describe how to the Drude peak arises within the memory function formalism in Section 3. 
Section 4 describes the holographic computation of cr(a;), and a proof of equivalence with the results of 
Section 3. 


2. Hydrodynamics and the Drude Peak 


Drude peaks generically emerge from quantum field theories described by hydrodynamics with a small 
amount of momentum relaxation. Unlike the historical model of Drude physics (electrons scattering off 
of a lattice or impurities), quasiparticles need not exist to observe (1). Consider the long-time equation 
for momentum relaxation in a charged fluid, placed in an electric held in the x-direction; 


dtU^ + dxP 


IJx 

- +QE,, 


( 2 ) 


Here Ux is the x-momentum density, P is the pressure, r is the momentum relaxation time, and Q is the 
charge density (at rest) of the huid. As we are interested in zero-momentum transport, at frequency ui 
this equation simplihes to 


- - iu; ] Ux = QEx- 


( 3 ) 


To compute 


^dc — 


Jx 


( 4 ) 


^This follows directly from the fact that at finite density, one can generate an electrical current without any energy in 
a translationally invariant theory by simply boosting to a reference frame with a relative velocity to the rest frame of the 
metal. 
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we need to relate Jx to TT^. 

At this point, let us specialize further to a relativistic quantum field theory, and work in units where 
h = c = 1? All holographic models we will study will fall into this class of theories. In such a theory, we 
find that at small velocities Vx, the momentum density is given by 

TIj, = (e + P)vx, (5) 

where e is the energy density of the fluid at rest and P is the pressure, and 


Jx — QVx- 


( 6 ) 


Combining these equations we can relate Jx to Hx, and find a Drude conductivity of the form (1) with 


<^dc — 


Q^t 
e + P' 


(7) 


This is a special limiting case of more general hydrodynamic results found in [26] . It is valid when charge 
diffusion is negligible, and there is no magnetic field. 


3. Memory Functions 


Hydrodynamics does not give us an explicit expression for any thermodynamic quantities, or r. The 
thermodynamic functions such as e and P are “intrinsic” to the quantum field theory, and are approxi¬ 
mately independent to the precise mechanism of momentum relaxation (so long as is small). On the 
other hand, r is “extrinsic” and sensitive to the precise way in which momentum can relax. The memory 
function analysis provides us with a way of computing r, given that we know the microscopic method by 
which momentum can relax. 

More precisely, in a time-reversal symmetric theory in which momentum is the only almost conserved 
quantity, the memory matrix framework tells us that 


o-(w) 


Xjp 


Mpp{u}) — iujxpp' 


( 8 ) 


where xpp is the momentum-momentum susceptibility, xjp is the current-momentum susceptibility, and 
Mpp{uj) is the momentum-momentum component of the memory matrix. If there are other long-lived 
conserved quantities, then a matrix generalization of this equation applies [23]. For us, Mpp{oj) is a small 
quantity and may thus be computed perturbatively. In our relativistic field theories, one finds 


XPP = e + -P, (9a) 

Xjp = Q- (9b) 


To derive these results, we note that the conjugate thermodynamic variable to momentum Px is velocity 
Vx- The susceptibility XaP^ is equal to {a)/vx as Vx —>■ 0. (9) then follows from (5) and (6). As w —)■ 0, 
one finds [7] 


Mpp(0) = lim 

oj^O 


LV 


e + P 

T 


( 10 ) 


^Note that, from the perspective of condensed matter physics, “c” here is the effective emergent speed of light in the low 
energy Lorentz-invariant theory. 
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where Px is the time derivative of the total x-momentum. We will argue in the course of our holographic 
discussion that finite frequency corrections to the memory matrix will be extremely small in the regime 
we are interested in. r may thus be computed from this memory matrix component at w = 0, and a 
simple calculation verifies that (1) and (7) are recovered. 

Now, let us specialize further to a quantum field theory with a translationally invariant Hamiltonian 
Hq, perturbed to 

H = Ho- j d'^x /i(x)0(x), (11) 

where O is a Lorentz scalar operator in the original theory, and h{x) is a (real-valued) static, spatially 
varying field. The quantum operator 

Px = ^[Px,H] = -i y h{x)[Px,0{x)] = j d'^x h{x){dxO){x) (12) 

We then employ 

^(k)h(q)G'|^oa^o(k,q,w) (13) 

The perturbative parameter that makes Mpp small is ~ h?. At leading order in h, is the Green’s 
function of a translationally invariant theory and so the integral identically vanishes when k + q ^ 0. 

If we have a “lattice” or periodic potential where 


h{x) = ho cos(ko • x), 


(14) 


then we obtain 


G 


R 




(15) 


PxPx ' 2 

In the case of disorder, the function h(x) will be random. It is common to take h{x) to be a zero-mean 
Gaussian random function with mean and variance given by (E[- • • ] denotes disorder averaging) 

E[/i(x)] = 0, (16a) 

E[h(x)h(y)] = e^6(x - y). (16b) 

We obtain [24] 


G 


R 

PxPa 


(cj) = J d'^kd''qE[/i(k)h(q)]A:2Ggc,(k,q,a;) = J d'’*kd'’'q £2/0^6(k-h q)Ggc,(k,q,a;) 

= [ d^e^klG^oiKu^)- (17) 


We approximate ^ by its average as fluctuations are suppressed in the large volume limit [10] . More 

-^x 

generally, if we assume that h{x) and its derivatives are non-vanishing (almost) everywhere in space, we 
will find 

(18) 


p (cu) “ = ” 

p xP X 


d^\^\h{\^)\^klGlo%u:). 


We have belabored this rather trivial discussion to emphasize the following - in the last equation, the 
quotes around the equals sign arise because technically, we have neglected a 6 function enforcing momen¬ 
tum conservation in j/i(k)j^.^ However, writing the equation in the form (18) is convenient: we will see a 
similar equation with a “missing” 5 function arise holographically. Thus 


e + P 


T 


^ f d-k X l.m 

J uj->-0 (jO 


(19) 


®For example, in the disorder case, using h(—k) = /i(k) and E[h(k)/i(q)] = e^5(k + q), (but without the 5 function) in 
(18), we obtain (17). Crudely, one should think of dividing |h(k)p by 5(0) in (18). 
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4. Holography 


Now let us turn to holography. We assume that d > 1. We consider solutions of the Einstein-Maxwell- 
dilaton (EMD) system with action 


s = y (di - ndutf - 


Z{t) 


FrsF^^^ , 


( 20 ) 


where gMN is the bulk metric (dual to the stress tensor of the boundary theory), Am is a U(l) gauge field 
(dual to the electric current of the boundary theory), and is a dilaton field. These holographic models 
are known to give rise to rich families of boundary theories; we focus on metallic phases in this paper. 

Let us now, without proof, state some results about static, isotropic geometries that solve the equations 
of motion associated with this EMD action. We use the conventions and results of [10] in what follows. 
These solutions have metric 


ds2 


— f^dr2 

r2 


a{r)b{r)dt‘^ + dx^ 


( 21 ) 


where 6(r) plays the role of an “emblackening factor”: in particular, near the black hole horizon (of planar 
topology), located at finite r = rh, we find 


b{r) « 4'7Tr(rh — r), 


( 22 ) 


where T is the Hawking temperature of the black hole, which also equals the dual field theory’s temper¬ 
ature; a is finite near the horizon. Near the boundary (r = 0), a(0) = 6(0) = 1; <P{r = 0) = 0, Z(0) = 1 
and E(0) = —d{d + 1); the asymptotic geometry is that of AdS. The UV of the continuum field theory is 
thus approximately conformal. The profile of the gauge field is 

A = p{r)dt (23) 


where p{rh) = 0, and the asymptotic behavior near the boundary is 


p{r) 




(24) 


/i is the chemical potential associated with the conserved charge. One finds that the object C, defined as 


(abY 2 Qk? 

ard i^d 


(25) 


is independent of r, where here and henceforth, primes denote r-derivatives. There is, of course, much 
more that can be said about these geometries, but this is all we will need for the present paper. Although 
we have assumed (for technical ease) that the UV geometry is AdS, the remainder of the geometry may 
be completely arbitrary, so long as it may be constructed as a solution of EMD theory. 

Now, we will add to this geometry a fourth field: a neutral scalar ■i/’, dual to the operator O sourced 
by the translation symmetry breaking field /i(x). The action of ip is 


= j {{dMpf + , 


(26) 


If the operator O has (UV) dimension A> (d-|- l)/2 (we choose this so that the operator O is described 
by standard quantization in holography), then the near-boundary asymptotic expansion of O is 


p{r) = r' 


4-1-1 




X 


Ld/2 


+ --- + B 


LV2 


+ 


(27) 
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and -B(O) = A[A — d — l)/Lp‘. In general, B should be chosen non-trivially in order to obtain results which 
match boundary theory expectations [9]. The boundary conditions which imply that the Hamiltonian of 
the boundary theory is given by (11) are imposed by setting 

=/i(x). (28) 

This will induce a small backreaction on the original EMD system, as the bulk stress tensor of V’ is a source 
in Einstein’s equations. However, this source is 0{h?), and following [8, 9] we may treat the geometry at 
0(/i) (thus the background metric is unperturbed) in order to compute the conductivity at leading order 
in h (as will be clear in the derivation shortly). We may thus take 

= j d\h(k)'iljo{k,ry^-^, (29) 


with '0o(k, r —)• 0) ~ j_^-d/2^d+i-A ^ obeying the linear equation of motion associated with the 

action (26), subject to a regularity condition at r = rh, and the boundary condition (28). A holographic 
calculation can explicitly determine the scale of h at which this approximation breaks down [9]. 

The expectation value of 0{x) in our state at finite density and temperature, and in the background 
field h, is given by [30] 

{0{x)) = {2A- (30) 
If the asymptotic expansion of near r = 0 is 

= + + , (31) 


then 


a{uj) 




We will use these facts to compute Green’s functions and conductivities in this section. 


(32) 


4.1. Drude Peak 


To compute the conductivity, we need to compute the response of our background EMD-scalar solution 
to a perturbation 5AxQ~^'^^ (at zero momentum) - equivalent in the boundary theory to imposing a small 
electric field. The solution need only be found within linear response theory, as is standard, but the 
computation below will proceed a bit different than that in [8, 9], as we must explicitly consider finite u) 
effects. 5Ax can only consider spin 1 perturbations under the spatial SO(d) isometry of the metric. These 
perturbations are (in axial gauge 5grx = 0): 


5Ax, Sgtx = j2^gtx, dxdii, d'^dxd'ip, ... 


By computing the full linear response problem (we discuss the remainder of the necessary boundary 
conditions later), we can find the near boundary asymptotic behavior of 6Ax{r) and compute cr(w) from 
(32). 

Linearizing the EMD system, one finds the equations of motion: 





(33a) 

(33b) 
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fcxC^V>o(k,r)^ 

brd 


(5V’(k, r) 
V’o(k,r) 


(33c) 


^9tx 


(^V’o(k,r)2 





+ ^V'o(k,r)5V'(k,r-), 


where we have defined 




b 

r'^uj 


d'^k k^\h{\^)\^1po{]^,rf 


( Sip{uj,k,r) 
V ^o(k,r) 


/ 


(34) 


In these equations, 5Ax, 5gtx and 6Vx are Fourier components at zero momentum and finite (but small) u; 
they are functions of r only. The equations of (33) are in order: the rx-component of Einstein’s equations, 
the x-component of Maxwell’s equations, and the ijj wave equation at momentum k; all of these equations 
have been integrated over all of space. Up to the O(cu^) terms in the scalar EOM, the equations close to 
a finite set of equations for 6Ax, 5Vx, and 5gtx- Since we want to work at finite frequencies w <C T, but 
UJT possibly S> 1, it requires some care to argue that we need only consider these three perturbations. 

We wrote coefficients of the form f d‘’*k |/i(k)p-i/’o(k, r)^ above. This is analogous to (18) - there is a 
“missing” 5 function. The reason this arises is that we have integrated over all of space in order to isolate 
the zero momentum mode SVx - this induces an “inhnity” factor analogous to 5(0) in any term which 
does not vary in space: for example, the other two terms in (33a). Properly regulating this inhnity, we 
recover the equations above - provided that |/i(k)p is understood to be “missing” a factor 5(0), as in 
Section 3. 

Our strategy is to show that the linearized huctuations need only be computed to 0(a;) in order to 
compute the linear response problem to all orders in cur. This will imply that we can focus on the w —?• 0 
limit of the equations of motion, where we need only solve a linear response problem in three variables, 
instead of an inhnite number. 

Let us begin by studying the linearized equations of motion when cj = 0 exactly - we will not worry 
about computing the conductivity just yet. In this case, the equations of motion reduce exactly to 


- Q'^Ax L<^5Vx, 



5V' 


bgtx 


1 

b^-d 




(35a) 

(35b) 

(35c) 


These equations are identical to those of massive gravity [13] in the limit a; —?■ 0, with the object in 
brackets in (35c) playing the role of the graviton mass. It is important to remember that it is 5Vx which 
stays 0(1) as w —)• 0, and not 5ijj{uj,k,r), which is 0(a;). Let us exactly find - at leading order in h - 
all of the solutions of these equations; counting derivatives we find there must be 4 linearly independent 
solutions. The hrst is spotted by inspection (we will not normalize any of these linearly independent 
solutions with the correct dimensions): 

5Ax = l, S'Px = -^, 6gtx = 0. (36) 

In fact, this is the only perturbation that, at leading order in h, couples to the '0 sector. All modes will 
couple to this sector, but only at 0{h?), which will prove to be subleading in our computation of the 
conductivity. The remaining three modes are, at leading order, modes of the translationally invariant 
theory with “ip = 0, and may also be written down exactly. Eirst there is a “diffeomorphism” mode 


5gtx = 1 , SAx = 0 , 5Vx = 0 . 


(37) 
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The next mode may be found by performing a Galilean boost to the static solution: 


6A^ =p + p^, 6gtx = l-ab, 5Vx = 0. 


( 38 ) 


where ph is a constant, and we have conveniently chosen that 5gtx{r = 0) = 0. We may fix the value of 
Ph by requiring that (35a) be satisfied at r = rh, and we find, using (22) and (25) (along with the fact 
that C = 0 is hxed by the linearized equations of motion): 


_ 2-71_ Ts 

K?r'^Q Q ’ 


(39) 


where s is the entropy density (we have converted between the area of the black hole horizon and the 
entropy density of the dual field theory using the Bekenstein-Hawking formula). We may use the reduction 
of order method^ [31] to compute the final solution, by using the fact that both (39), and the mode we 
are looking for, solve (at leading order) the differential equation 


and we find 


2e^K?Q^ar^ 

JJ2d-2 


6 Ax = br^-'^Z6A' 


(40) 


b{s)Z{s){p{,) + SV, = 0. (41) 

0 0 

Note that, at leading order in h, this mode has a logarithmic divergence at r = in bAx alone. 

As (Tdc must be hnite, and we must have bgtx{r = 0) = 0 (we are not sourcing any temperature 
gradients), at w = 0 the solution to the linearized equations of motion obeying all boundary conditions 
turns out to be simply 

bAx = bAl, 6Vx = ^bAl 5gtx = 0. (42) 

We now wish to include corrections when u <C T, p, but work to all orders in lot. More precisely, we 
set all subleading corrections in uj/T or co/p to vanish, and then to all orders in lot. First, we need to 
discuss the boundary conditions at the black hole horizon. Both bVx and 5Ax obey infalling boundary 
conditions: 

M., bVx ~ (rh - (43) 

These boundary conditions arise from including the 0(0;^) terms in the equations of motion, and demand¬ 
ing that perturbations fall into the black hole horizon (in the boundary theory: energy is dissipated). They 
are, at first glance, non-perturbative in ui, which is frustrating for our purposes. However, things are not 
so bad. Let us fix rh — r to be arbitrarily small, but send a; —?■ 0. We may then Taylor expand 


= 1 + ^ log ^[1^ + . . . (44) 

d-TTi rh — r 

In a neighborhood of this fixed r, the Taylor expansion in to must be a solution of the equations of 
motion, but the w-dependence of the equations of motions themselves only comes in at 0(a;^). We 
therefore conclude that the 0(a;) coefficient must be a solution to the equations of motion evaluated at 
a; = 0. As we are going to argue later that it is sufficient to only compute 6Ax, bgtx and bVx to 0(t<;), 
this is a great simplification - as far as our computation is concerned, (35) is exact. Imposing infalling 

^Suppose we have a differential equation fay + (fiy')' = 0, and yo{r) is an exact solution of this differential equation. 
One can show that a linearly independent solution is yi{r) = yo{r) dsfi{s)~^yo{s)~^. 
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boundary conditions becomes equivalent to imposing the boundary conditions (up to 0(a;) terms finite 
at the horizon): 


6A^{r rh) = 


1 + 


ltd 


log 


rh 


Vh-r 


+ ... 


5V.{r 




1 + 


lUJ 


log 


rh 


rh — r 


+ 


(45) 


We impose no specific boundary condition on Sgtx at r = rh, other than no logarithmic divergences. 

Let us now go ahead and compute the 0(a;) corrections to the linearized modes, at leading order in 
h. Imposing infalling boundary conditions on 6Vx, (35c) implies 


iw Q 


5^0 


1 

rh-r 


^9tx 

47rTr^(rh — r) 


d‘^k|fi(k)|^A:^V’o(k,rh)^. 


(46) 


Evidently, this equation is only consistent if 6gtx{r = rh) ~ is a finite number. Since our boundary 

conditions are that 5gtx{r = 0) = 0, we conclude that 5gtx is not a constant at leading order. (35b) then 
implies that 5Ax must have a component at 0(/i“^) as well. It is consistent in every equation, at 0{uj), to 
take 5Vx to be 0{hP), so long as the 0{uh ^) coefficients in 6Ax and Sgtx obey the equations of motion 
associated with the translationally invariant black hole. In fact, at 0(a;/i“^), we hnd that the boost mode 
(39) is the only mode consistent with all of our boundary conditions: (37) is ruled out by 5gtx{r = 0) = 0, 
and (41) by the fact that the logarithmic divergence in 5Ax occurs at 0{hP). Using (46) we in fact find 
that the coefficient of the Galilean boost mode in 5Ax is: 

M,, = MS - ^(P + Ph)6Al + O (e°) , (47) 

where we have defined 

^ /d'^k|/i(k)p/cSV’o(k,rh)^. (48) 

r r^ J 

In the next subsection, we show that this r is equivalent to that defined via the memory matrix; for now 
this is simply a definition of r. Note that the logarithmic divergence in 5Ax is not included, as it is 0{hP). 

Now, let us explain why we do not have to worry about any further terms, in so far as computing 
the conductivity at leading order in our limit. Firstly, let us consider the term that arises from the 
logarithmic divergence at 0{hP) in SAx- This will indeed lead to subleading corrections to the Galilean 
mode described above; however, these corrections, by dimensional analysis, will scale as {oj/T)F{T/g)5A^ 
for some scaling function F, and are subleading in our limit. A similar argument holds for the 
term which arises from the fact that when ujt S> 1, the Galilean mode dominates the constant contribution 
5A^: the correction due to infalling boundary conditions on this term ~ w^r/T, and this can be neglected 
in our limit. Similar arguments will hold for all further corrections at higher orders in ca. In a nutshell, the 
linearized equations contain perturbations at most ~ which implies that only the terms need 

to be included at leading order. The fact that the boost mode is the only important 0(a;) contribution 
to 5Ax is suggestive of the fact that this is “hydrodynamic” transport. 

To compute <j{oj) we simply employ (24) to find the near boundary behavior of p{r). We thus find 


M(0)(a;) = M' 


iwQr / Ts\ 




1 - 


iujT^pQ + Ts) 
e + P 


= - ^^r), 






(49a) 

(49b) 


where we have used the thermodynamic identity 


e + P = pQ + Ts 


(50) 
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in the last step of (49a). Recall that /5 is related to cr(w) by (32). It is straightforward from 
here to obtain (1), (j{uj) = crdc/(l — iwr), and the expression (7) for a^c in tbe limit where u/T —)• 0 hrst, 
and ujT is held finite. 

That Drude physics arises is not particularly surprising for a few reasons. Intuitively, we expect 
via the fluid-gravity correspondence [32] that a holographic system, weakly perturbed by translational 
symmetry breaking, behaves in the same way as a fluid with momentum relaxation would: see also [33]. 
Additionally, it was shown via matched asymptotic expansions in [12] that certain massive gravity theories 
contain Drude peaks in the limit of weak graviton mass, at T = 0. Due to the equivalence between the 
holographic set-ups above and massive gravity as w —>■ 0 [8] (with the caveat that the graviton mass 
becomes a function of r), this result makes sense. It is pleasing nonetheless to see the Drude peak emerge 
from a direct calculation for a wide range of holographic theories at finite T. 


4.2. Equivalence with the Memory Function Approach 


Our only remaining task is to show that the r defined above is equivalent to the r defined in (19), defined 
in terms of the leading order imaginary behavior of G^q{uj —)• 0). By the holographic dictionary for 
retarded Green’s functions, we impose infalling boundary conditions at the black hole horizon r = rt, as 
before. Near the black hole horizon: 


'(/’(k, r, w —)■ 0) •i/;(k, r, a; = 0) 


1 + 


lUJ 

4)^ 


log 


r\y 


rh - r 


(51) 


As before, at 0(a;), 'll: solves its cu = 0 equation of motion. Also note 'ip{'k,r,u} = 0) = with i/jq 

the finite scalar profile we defined previously. We wish to “propagate” this small imaginary piece from the 
horizon to the AdS boundary. This directly allows us to compute Im(G^(a;))/a; by using the AdS/CFT 
dictionary at r = 0. We compute the linearly independent solution to equation of motion with the 
reduction of order technique: 


V’o(k,r) 


ds 


5(s)V’o(k, s)2 


(52) 


We assume that 'i/jq has the same asymptotics as before: V’o(k, r —)• 0) ~ ]^-d/ 2 ^d+i-A^ have normalized 
■01 conveniently. The integral as written above is convergent so long as A > (d-|-l)/2, and we have assumed 
this property previously. In fact, this solution also has the r —?■ 0 asymptotics we wish: as 6 ~ 1 near the 
AdS boundary, we hnd that 


V^i(k,r)«L-^/2- 


(r —)• 0). 


(53) 


2Z\ - d - 1 

Near the black hole horizon, 'ipi is divergent - this divergence implies that the imaginary contribution 
to ^|J is (for the purposes of our calculation) proportional to 'ipi. As this divergence is associated with 
near-horizon physics, its coefficient is completely independent of the UV or intermediate scales in the 
geometry: 


'0i(k,r) 


log- 




4'iTrV’o(k, rh)L^ Th — r 
Comparing (51) and (54) we conclude that to order w, 

/o \ I n ^ ■ -^^%(k,rh)2 

V’(k, r, uj) ^ ipoiK r) + iw-^- 


-|-hnite, (r —)• rh). 


V’i(k,r)-FO(w2), 


The function is real. Therefore using (30), it is easy to conclude that 

Im(G ggi(k,u;)) _ L%(k,rh)2 
OJ 


lim 


(54) 


(55) 


(56) 
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Comparing (19), (48) and (56), we see that the conductivity of this metal, computed directly via holog¬ 
raphy, is identical to the result computed with the memory function formalism, without using any holog¬ 
raphy. These results are valid in the limit when are vanishingly small compared to T and /r, the 

regime of validity where the approximation (8) holds. 


4.3. Seebeck Coefficient 

Let us briefly discuss thermoelectric transport coefficients [34, 35, 36, 37, 38]. One such coefficient is 
immediately computable given our discussion of the conductivity: the Seebeck coefficient a, defined by 


Qx — OiTExi 

where Qx is the heat flow density in the x direction. Using 

qx = {T^n-f^{Jl, 


(57) 


(58) 


and reading off the expectation value of (T**) in the linearized modes computed in this section, it is 
straightforward to recover 


a{uj) 


sQt 1 

e + Pl-\u)T' 


(59) 


in agreement with hydrodynamics [26] and the memory function formalism (use that the heat-momentum 
susceptibility xqp = sT). We also expect that a calculation of the thermal conductance it, defined by 


Qx — 

would find, in agreement with hydrodynamics and the memory function formalism. 




s^Tt 1 


(60) 


(61) 


e + Pl-\0JT 

Computations of R{uj) are a bit more involved [34, 35, 36, 37, 38] and we will not pursue them further. 


5. Conclusion 


In this work, we have demonstrated an exact correspondence between the memory function approach 
and a holographic approach to transport in a strongly correlated “strange” metallic phase of matter 
without quasiparticles, in the limit of finite density and slow momentum relaxation. This result crispens 
the qualitative scaling agreements noted in special cases earlier [8, 9]. We have also pointed out the 
emergence of a universal Drude peak in these holographic models. These results are to be expected on 
physical grounds, but it is nonetheless instructive to see them explicitly verified. 

The reason that we add a fourth scalar held V’) instead of say adding perturbations to the dilaton, is 
that the EMD background is only corrected at 0(/i^); it would receive corrections at 0{h) if the dilaton 
had a component which breaks translational symmetry. As in [8, 9], adding a new scalar simplihes the 
resulting computation of 17 ( 0 ;). It would be worthwhile to explicitly show that the correspondence with 
the memory matrix survives even if we break translational symmetry through one of the EMD helds: for 
example, using results of [6], it should be possible to do this for spatially-dependent chemical potentials. 

One of the advantages to a holographic computation is that we are formally not restricted to the 
limit of weak momentum relaxation. Approaches including massive gravity and Q-lattices, which produce 
translation invariant geometries, allow for analytic control in this limit. Perhaps techniques similar to 
ours can remain valid even in this limit of strong momentum dissipation. While it is now clear that such 
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approaches mimic more conventional “lattices” or disorder when momentum relaxation is weak, it is less 
clear whether this analogy remains when momentum relaxation is strong. It is possible that more exotic 
phases of holographic matter, analogous to a quantum glass or many-body localized phase [39], will arise 
when translational symmetry is explicitly broken strongly. Alternately, efficient momentum relaxation 
without “glassy” physics may be responsible for scaling properties of cuprate strange metals [40] ; perhaps 
massive gravity/Q-lattice approaches are powerful tools for these problems. More work in this direction 
is warranted. 
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